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Abstract 


Classical molecular dynamics (MD) simulations employing Brenner po- 
tential for intra-nanotube interactions and van der Waals forces for polymer- 
nanotube interface have been used to investigate thermal expansion and diffu- 
sion characteristics of carbon nanotube-polyethylene composites. Addition of 
carbon nanotubes to polymer matrix is found to significantly increase the glass 
transition temperature Tg, and thermal expansion and diffusion coefficients 
in the composite above Tg. The increase has been attributed to the temper- 
ature dependent increase of the excluded volume for the polymer chains, and 
the findings could have implications in the composite processing, coating and 
painting applications. 
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Carbon nanotubes (CNTs) are nano-scaled materials that are found to have unusual 
mechanical, electronic and electromechanical properties. 1 The strong in-plane graphitic C- 
C bonds make them exceptionally strong and stiff against axial strains and very flexible 
against non-axial strains. Additionally, CNTs can be metallic or semiconducting, deter- 
mined by their atomic configurations such as chirality and diameter. Many applications of 
CNTs, such as in nano-scale molecular electronics, sensing and actuating devices or as rein- 
forcing additive fibers in functional composite materials have been proposed. Several recent 
experiments on the preparation and mechanical characterization of CNT-polymer compos- 
ites have also appeared. 2-4 These measurements suggest modest enhancement in strength 
characteristics of CNT-embedded matrix as compared to the bare polymer matrix. 4 Pre- 
liminary experiments and simulation studies on the thermal properties of carbon nanotubes 
show very high thermal conductivity. 5 It is expected, therefore, that nanotube reinforcement 
in polymeric materials may also significantly change the thermal and thermo-mechanical 
properties. 

Using classical molecular dynamics (MD) simulation, we investigate thermo-structural 
behavior of CNT reinforced polyethylene (PE) composite. The structural and diffusion 
behavior of the composite has been studied in the temperature range below and above 
glass-transition temperature Tg. We report significant enhancement of thermal expansion 
and diffusion coefficients of the composite material above Tg. 

The interactions of the atoms within the nanotube are described by Tersoff-Brenner 
many-body potential, 9,10 and the nanotube-polymer interface is described via non-bonding 
van der Waals interactions (VDW) of truncated Lennard-Jones (LJ) 6-12 type. The dy- 
namics of polymer chains is described by a united-atom PE model (DLPOLY package 7 ), in 
which each CH2 or CH3 groups are considered as single interaction units. The interactions 
within a polymer chain have two components 8 , valence angle potential: 

<£(#) = O.5&0 (cos 9 — cos # 0 ) 2 
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where k e = 520 kJmol~\ 6 0 = 112.813°, and torsional potential: 

$(&)/ — Co + C\ coso! + C 2 Cos 2 q + C 3 COS 3 q: 

where C 0 = 8832, C x = 18087, C 2 = 4880, and C 3 = -31800. The C-C bonds in a polymer 
chain can be kept rigid to have large time steps. However, we allow full dynamics of C- 
C bonds, to be consistent with the C-C bonds within CNTs, and use harmonic potential 
0.5 kb{l — lo) 2 for the C-C bonds within a chain, where k b = 346K J mol -1 A . n Truncated 6- 
12 LJ potentials are applied to pairs of units parted by more than three units. The time-step 
0.5fs is used in MD simulations. 

The conformation of long-chained polymer molecule is important in describing the dy- 
namics and properties of bulk polymers. 12 Studies have shown that MD is very slow to relax 
a polymer chain to equilibrium conformation. This difficulty can be avoided by using Monte 
Carlo (MC) simulations. In our study MC simulations were run on each chain at 300K for 
up to two millions of steps, such that the end-to-end distances of polymer chains show the 
square-root dependence on the chain length, in agreement with the predictions of Flory s 
theory. 13 MC relaxed chains were then used in MD equilibration of bulk configurations with 
periodic boundary condition. The initial overlapping of chains is removed by gradually 
turning on of the VDW interactions. 

We first investigate a bulk PE sample of 80 short chains, with 10 units (of about 10 
A long). For nanotube reinforced composite samples a 20 A long capped (10,0) CNT is 
embedded in the polymer matrix (Figure 1), with about 8% volume ratio. The Evans 
NVT 14 MD was run for 100 ps followed by lOOps Berendsen NPT 15 (P = lbar, T= 300K) 
MD to help samples reach their initial equilibrium configurations. The samples are then 
gradually cooled down to 10K at a rate of IK/ps. Each system is further equilibrated for 
100 ps at chosen temperature data points (at P = lbar) to reduce possible fluctuation. 

The temperature-depended densities of both the bulk polymers and polymer-CNT com- 
posites are shown in Figure 2. The location of the discontinuity in the slope of the density vs 
temperature plot implies the on-set of the glass transition temperature (Tg). Depending on 
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the molecular weights and degree of cross-linking, experiments have reported Tg for bulk PE 
varying from 180K to 250K. 16 In our sample, the computed Tg of about 150K is in agreement 
with previous MD simulations on similar size samples. 11 The nanotube-polymer composite, 
as shown in Fig. 2, tends to have a higher Tg of 170K, as the CNT in the matrix tends 
to slow the motions of the surrounding molecules below Tg. The cross-linking of polymer 
matrix with embedded CNTs can also further reduce the motions of polymer molecules and 
increase the Tg. Recent experiments on epoxy with 1% weight CNT composite indicates 
about 10 °C increase of Tg. 17 

The decrease in the density of the bulk polymer and the composite with increased tem- 
perature indicates thermal expansion of the materials. The slope of the curve at a given 
temperature gives thermal expansion coefficient at that temperature. It is noted that above 
Tg (170K), the thermal expansion increases more rapidly in the composite compared to that 
in bulk polymer. The volume thermal expansion coefficient of the composite is found to be 
ToxHT 4 /^ 1 below Tg, about 18% larger than the 3.8xlO- 4 Ff- 1 of the polymer bulk 18 in 
the same temperature range. Above Tg, the thermal expansion coefficient of the composite 
increases as much as 40% compared to the bulk polymer within the same temperature range. 

The simulations on the larger systems verified this increase of thermal expansion coeffi- 
cient of the composites. Samples with 50 PE chains of length of about 100A, with a 200A 
long (10,0) CNT embedded in the composite, were also studied. Figure 3 shows a similar 
increase of the thermal expansion in the composite above Tg, which is 16.1xl0“ 4 A' _1 , a 
increase of 142% compared to that of a bulk polymer. The Tg for the matrix is found to be 
around 300K, which is higher than Tg for small samples as discussed above. This is due to 
the dependence of the Tg on molecular weight of polymer chains. Previous MD simulation 
on long-chained PE systems show Tg to be in the range of 230K to 350K, 8 - 19 which is similar 
to our computed Tg for long chain systems. 

The origin of higher thermal expansion coefficient in the composite is explained as fol- 
lows. The CNT embedded within a polymer matrix has a fixed volume and excludes the 
occupancy of polymer chains. The enhanced thermal expansion of the composite is at- 
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tributed to the increase in the excluded volume as a function of temperature. To verify this, 
the MD simulations of the smaller-sized composite sample at 300K were repeated under two 
constraints. First, the phonon vibrations within a nanotube were frozen by increasing the 
stiffness of C-C bonds and allowing only the motion of the CNT as a whole. Second, in 
addition to the first, the nanotube as a whole was frozen as well. Results under the second 
constraint, for a completely frozen CNT, show that the density of the composite reverts 
back. The density of the composite under the first constraint is intermediate between the 
completely frozen and the completely free case. These results (Fig. 4) suggest that both 
phonon modes and Brownian motions of the CNT contribute to the excluded volume. The 
increase of the thermal expansion coefficient suggests that the linear mixing rule for macro- 
scopic fibers composites is not valid for CNT composites, because the linear mixing rule does 
not take into account of effects of thermal motions of CNTs on their surrounding polymer 
molecules in a composite. 

Lastly, in Fig. 5, we show the diffusion coefficients of the C atoms of the polymer 
matrix for the bulk polymer and composite cases. It is clear that in the composite case, the 
polymer chains are more diffusive at temperatures above Tg. Further more, the increase in 
the diffusion coefficient parallel to the CNT axis is about 30 percent larger than the increase 
perpendicular to the tube axis. This correlates well with the higher thermal expansion 
coefficient (above Tg) data described above, and shows that the increased phonon vibrations 
and Brownian motion of the CNT are not only coupled efficiently to the motions of the 
polymer chains but also the coupling is anisotropic. This means that polymer chains parallel 
to the tube axis will flow or diffuse better at higher temperatures. 

The importance of these findings have implications in the processing and applications of 
CNT reinforced polymer-composite materials. The higher diffusion coefficient of the matrix 
atoms above Tg will allow increased mobility of composite material during processing steps 
such as continuous spinning, weaving and extrusion of the material for fabrication purposes, 
and also possibly the smoother flow through nozzles for painting or coating applications. 
The increased thermal expansion and diffusion coefficients, above Tg, may make it possible 
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to significantly increase the thermal conductivity of the composite at high temperatures. 

In summary, we have used MD simulations to study temperature dependent structural 
behavior of polymer-CNT composites. The results show that the thermal expansion coef- 
ficient of the composite and diffusivity of polymer molecules increases significantly above 
glass transition temperature of the composite. 
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FIGURES 



FIG. 1. The MD simulation box for a composite system with 10 A long capped (10,0) CNTs 
embedded in polyethylene matrix 
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FIG. 2. The density as function of temperature for short-chained PE and short (10,0) CNT 
composite system. Tm is melting temperatures marked according to empirical formula Tm/Tg' 
1.5. Systems are in a rubber-like state between Tg and Tm 
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FIG. 4. The change on density of CNT composite with phonon modes and Brownian motion 


frozen. Both thermal motions contribute to the exclude volume of embedded CNTs, which lead to 
a increased thermal expansions. 
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FIG. 5. The diffusion of polymer chains. With the presence of CNTs, the diffusion coefficient 
is increased at high temperatures ( Tg) and the components parallel to CNT axis has a larger 
contribution compared with the perpendicular one. 
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